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Abstract 

Using the dynamical mean-field theory, we calculate the effective electron mass in the Hubbard model on a 
semi-infinite lattice. At the surface the effective mass is strongly enhanced. Near half-filling this gives rise to a 
correlation- driven one-electron surface mode. 



Important characteristics of interacting elec- 
trons on a lattice can be studied within the Hub- 
bard model. If symmetry-broken phases are ig- 
nored, the electron system is generally expected 
to form a Fermi- liquid state in dimensions D > 1. 
The effective mass m* which determines the quasi- 
particle dispersion close to the Fermi energy, can 
be substantially enhanced for strong interaction 
U. A further enhancement may be possible at the 
surface of a D = 3 lattice since here the reduced 
coordination number z s of the surface sites results 
in a reduced variance A oc z s of the free local 
density of states. This implies the effective inter- 
action U I VA to be larger and thereby a tendency 
to strengthen correlation effects at the surface. 
Using the dynamical mean-field theory (DMFT) 
jjj , we investigate the strongly correlated Hubbard 
model with uniform nearest-neighbor hopping t on 
a semi-infinite sc(100) lattice at zero temperature. 



Close to half-filling, the surface effective mass is 
found to be strongly enhanced compared with the 
bulk value. 

We have generalized DMFT to film geome- 
tries A film consisting of d = 15 lay- 
ers with the normal along the [100] direction 
turns out to be sufficient to simulate the actual 
sc(100) surface, i. e. bulk properties are recov- 
ered at the film center. The lattice problem 
is mapped onto a set of d single-impurity An- 
derson models (SIAM) which are indirectly cou- 
pled via the self-consistency condition of DMFT: 
G { Q a) (E) = {G aa {E)- x +T, a {E))~ l . The itera- 
tive solution starts with the layer-dependent lo- 
cal self-energy S Q (a = 1, ...,d) which determines 
the layer-dependent on-site elements of the one- 
electron lattice Green function G aa via the (lat- 
tice) Dyson equation. The DMFT self-consistency 
equation then fixes the free impurity Green func- 
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Fig. 1. m* for the top- and sub-surface layer and the 
bulk as a function of the filling n. U = 24|f| = 2W 
(W with of the free bulk density of states). 



tion Gq"' of the a-th SIAM. Following the exact- 
diagonalization approach H , we take a finite num- 
ber of n s sites in the impurity model. The param- 
eters of (the a-th) SIAM are then obtained by 
fitting Gq*^ . Finally, the exact numerical solution 
of the a-th SIAM yields a new estimate for S Q 
which is required for the next cycle. From the 
self-consistent solution the effective mass is calcu- 
lated as to* = 1 - dE a (0)/dE. 

The figure shows the bulk mass m£ to increase 
with increasing filling n. For U = 2W the system 
is a Mott-Hubbard insulator at half-filling (n = 1). 
Consequently, m* diverges for n i— ► 1. Up to n w 
0.98 all layer-dependent masses to* (a — 2,3, ...) 
are almost equal except for the top-layer mass m* 
(a = 1). For n i— > 1, to* is considerably enhanced, 
e. g. m*Jm* b = 2.4 for n = 0.99 @. 

There is an interesting consequence of this find- 
ing: Close to the Fermi energy where damp- 
ing effects are unimportant, the Green function 
G a p(k,E) can be obtained from the low-energy 
expansion of the self-energy: T, a (E) = a a + 
(1 — rn* a )E + ■ ■■. For each wave vector k of the 
two-dimensional surface Brillouin zone (SBZ), the 
poles of G a p (k, E) yield the quasi-particle energies 



f] r (k) (r = l,...,d). The r) r (k) are the eigenval- 
ues of the renormalized hopping matrix T a/ g(k) = 
(m a mf])~ 1 / 2 (e a f3('k)+5 a f3(a a — fi)) where the non- 
zero elements of the free hopping matrix are given 
by eii(k) = e QQ (k) = 2t(cosfc x + cosfcj,) and 
ej_(k) = e Qa ±i(k) = \t\. For the semi- infinite 
system (d i— > oo) the eigenvalues form a contin- 
uum of (bulk) excitations at a given k. A sur- 
face excitation can split off the bulk continuum if 
m*/ml is sufficiently large. Assuming m* = to*. 
for a / 1 and to* =1 = m* (cf. Fig. 1), a simple 
analysis of T a g(k) yields the following analytical 
criterion for the existence of a surface mode in 
the low-energy (coherent) part of the one-electron 
spectrum (r 2 = to^/to* < 1): 
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The evaluation for the present case shows that 
indeed a surface mode can split off at the (0, 0) 
and the (tt, it) points in the SBZ for fillings 
n > 0.94 and n > 0.90, respectively. The surface 
modification of a a — S a (0) (second term on the 
r.h.s) turns out to be unimportant for the effect. 
The excitation found is thus a correlation-driven 
surface excitation which is caused by the strong 
specific surface renormalization of the effective 
mass. 
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